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Abstract 

In the context of the interacting boson model with s, d and g bosons, the conditions 
for obtaining an intrinsic shape with octahedral symmetry are derived for a general 
Hamiltonian with up to two-body interactions. 
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1 Introduction 


The collective model of the atomic nucleus assumes a description in terms 
of the shape of a surface and its oscillations around that shape. In nuclei 
quadrupole deformations with multipolarity A = 2 parameterize the most 
important deviation from spherical shape. Their manifestation, either as a vi¬ 
bration around a spherical shape or through the mechanism of spontaneous 
symmetry breaking as a permanent deformation, is by now a widely accepted 
feature of nuclear structure [1]. Superimposed on quadrupole deformed shapes, 
octupole (A = 3) and hexadecapolc (A = 4) deformations are frequently con¬ 
sidered as well, usually in terms of vibrational oscillations although there is 
evidence for nuclei with a permanent octupole deformation [2], 

The accepted paradigm is that the spherical symmetry of the nucleus is re¬ 
duced to a lower one which in most cases corresponds to a shape with a single 
axis of symmetry. In technical terms the spherical SO(3) symmetry is spon¬ 
taneously broken and reduced to SO(2). Both SO(3) and SO(2) are examples 
of symmetries that depend on continuous variables, e.g., the three Euler an¬ 
gles or the angle of rotation around the axis of symmetry. Nothing prohibits 
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the further reduction of the symmetry of the intrinsic shape of the nucleus 
until a discrete (or point) symmetry with a finite number of invariance oper¬ 
ations remains or even until no symmetry remains. A well-known possibility 
is qnadrnpole deformation with three axes of different length. The intrinsic 
nuclear shape in that case has no continuous symmetry but is invariant under 
the eight transformations that form the discrete group D 2 h- (For a discussion 
of discrete groups see, e.g., Hamermesh [3].) The invariance of the nuclear 
shape under rotations and reflections—or the absence thereof in the case of 
deformation—determines the energy spectrum and leads to predictions that 
can be tested experimentally. 

In 1994 Li and Dudek [4] pointed out that intrinsic shapes with a higher-rank 
discrete symmetry can be obtained in the context of the collective model and 
that, specifically, an octupolc deformation with non-zero p = 2 component 
(and all other multipoles zero) exhibits the tetrahedral symmetry Tj. In sub¬ 
sequent studies [5,6,7,8] this observation was followed up systematically and 
it was shown how the tetrahedral, octahedral and icosahedral discrete sym¬ 
metries, Td, Oh and Ih, arise through combinations of deformations of specific 
multipolarity. 

For a fermionic quantum system, e.g., an odd-mass nucleus or a single-particle 
nuclear Hamiltonian, time-reversal invariance must be considered as well. 
Due to Kramers’ theorem, stationary eigenstates of such systems are at least 
twofold degenerate [9]. Inclusion of the operation of time reversal into the 
higher-rank discrete symmetries considered above leads to the enlarged ver¬ 
sions T°, 0° and 1° (see § 99 in Landau and Lifchitz [10] who denote these 
enlarged discrete groups with a prime). A prominent consequence of discrete 
symmetries of higher rank is the occurrence of more than twofold degenerate 
states in the spectrum of a single-particle nuclear Hamiltonian [4]. 

While the formal possibility of nuclear shapes with higher-rank discrete sym¬ 
metries is by now well established, the question remains whether such exotic 
deformations are realized in nuclei. Over the years this question has been stud¬ 
ied from a theoretical [11,12,13] and an experimental [14,15] point of view. The 
former studies have been consistently carried out in a mean-field approach 
usually supplemented with pairing correlations. Concerning the experimental 
studies, it is fair to say that up to now no conclusive evidence has been found 
that unambiguously establishes the existence of a nucleus with a higher-rank 
discrete symmetry. A notable example is the study of Jentschel et al. [16] who 
failed to find a vanishing quadrupole moment of a negative-parity band in 
156 Gd which should have been the ‘smoking gun’ of tetrahedral deformation. 

On the other hand, discrete symmetries are rather well established in light 
nuclei in connection with alpha-particle clustering which itself has a long his¬ 
tory in nuclear physics [17]. Algebraic models have been developed by Bijker 
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and Iachcllo [18,19] with the aim to describe the discrete symmetries D 3 h and 
Td associated with alpha-particles in a triangular or tetrahedral configuration. 
This approach has attracted renewed interest in view of novel experimental 
evidence recently found in 12 C and 16 0 [20,21], Discrete symmetries associ¬ 
ated with alpha-particle clustering, while an important and attractive field of 
activity, differ from those considered here where they arise in the context of 
the collective model of the nucleus. 


The aim of this series of papers is to analyze the question of the possible 
occurrence of higher-rank discrete symmetries in nuclei from a different the¬ 
oretical perspective. Since the work of Arima and Iachcllo [22] it is known 
that an alternative description of collective states in nuclei exists in terms 
of bosons in the context of the interacting boson model (IBM). Quadrupole 
collective states require s and d bosons, with angular momentum t = 0 and 
i = 2, respectively, and lead to the most elementary version of the model, 
the sd-IBM. Many refinements of this original version are possible [23] and 
already in the early papers on the IBM an / boson (£ = 3) is added to deal 
with negative-parity states with octupole collectivity [24,25,26]. Hexadecapole 
states, on the other hand, require the consideration of a g boson with i = 4. 


It will be shown that the higher-rank discrete symmetries, as encountered 
in mean-field approaches, can also be realized in the context of an algebraic 
model with the relevant degrees of freedom and that, for example, tetrahe¬ 
dral and octahedral symmetries can be studied in the context of the s/-IBM 
and sg-IBM, respectively. Due to the pervasiveness of quadrupole collectiv¬ 
ity in nuclei, the addition of a d boson will make these algebraic models less 
schematic, leading to the sdf -IBM and sdg- IBM, respectively. (In the former 
case it might even be necessary to consider the spdf-lBM with an additional 
negative-parity p boson with i = 1 [27].) 


This series starts with an investigation of octahedral symmetries in the frame¬ 
work of the sdg-IBM, adopting the model’s most general Hamiltonian with up 
to two-body interactions. The collective parameters of quadrupole and hex¬ 
adecapole deformation and their relation to octahedral symmetry are recalled 
in Sect. 2. In Sect. 3 the general Hamiltonian of the sdg-IBM is defined and its 
corresponding classical limit in the most general coherent state is derived. The 
catastrophe analysis of the resulting energy surface is carried out in Sect. 4 
with particular attention to the occurrence of minima with octahedral sym¬ 
metry. First conclusions from this analysis are drawn in Sect. 5. 
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2 Quadrupole and hexadecapole shapes 


Shapes with octahedral discrete symmetry occur in lowest order through a 
combination of hexadecapole deformations Y All (9, 0) with different //. To make 
the shape more realistic for nuclei, a quadrupole deformation should be added 
since that deformation is of lowest order in the geometric model of Bohr. With 
both these deformations the nuclear surface is parameterized in the following 
way: 


R(6, 0) = Ro(l + a 20 Y 2 o (0,0) + a 22 [Y 2 _ 2 (9, 0) + Y 2+2 (9, 0)] 

+ ®4C0 40 {01 0) + ®42[^4—2 1 0) + Y i+ 2 (9, 0)] 

+ O 44 [f 4—4 ( 9 , 0) + T 4+4 (M)]). (l) 

It is customary to define quadrupole-deformation variables through 


a 20 = 02 cos 72, 


a 2 2 


/T 
V 2 


02 sin 72, 


( 2 ) 


where 0 2 quantifies the quadrupole deviation from a sphere and y 2 the devia¬ 
tion from a quadrupole shape with axial symmetry. Similarly, a hexadecapole 
variable 04 is introduced which quantifies the deviation from a sphere. For the 
parameterization of hexadecapole asymmetric shapes two approaches have 
been adopted. In the first, a single 7 is introduced which parameterizes the 
deviation from axial symmetry at once for the quadrupole and hexadecapole 
degrees of freedom [28]. A symmetry argument then leads to a simplified pa¬ 
rameterization in terms of three variables 0 2 , 04 and 7 , 

°40 = y /^04 (5 cos 2 7 + 1 ), a 42 = y / || 0 . 4 sin 27 , a 44 = ^04 sin 2 7 . ( 3 ) 


In fact, there are three different such parameterizations [28] but Eq. (3) is the 
one that has been used up to now in the analysis of the sdg -IBM (see below). 
The parameter ranges are 0 < 0 2 < + 00 , — cx)< 04 <+oo and 0 < 7 < 7 t/ 3 , 
and the values 7 = 0 and 7 = 7 t /3 lead to a shape with axial symmetry. 

In a second approach, developed by Rohozinski and Sobiczewski [29], a sepa¬ 
rate, independent asymmetry parameter 74 is introduced for the hexadecapole 
deformation. Furthermore, to describe the full range of possible hexadecapole 
deformations, an additional variable S 4 is needed, which represents the convex¬ 
ity or concavity of the shape. This leads to five shape variables, two quadrupole 
and three hexadecapole ones which are related as follows to the original pa¬ 
rameterization ( 1 ): 
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040 = At {\Pk COS ^4 + \f^2 S ^ n COS 7 4 ) ’ 

042 = —\[\P4 sin 84 sin 74, 

044 = 04 {\pk COS $4 — \pk sin 84 COS 74^ , 


( 4 ) 


where the parameter ranges are now 0 < At < + 00 , 0 < 74 < 7 t /3 and 
0 < 84 < n. An axially symmetric shape occurs for 042 = 044 = 0 . This 
corresponds to 74 = 0 and $4 = arccos \J 7/12. 

A shape with octahedral symmetry implies a vanishing quadrupolc deforma¬ 
tion, a 2 o = o 2 2 = 0 , and can be realized in lowest order with a hexadecapole 
deformation that satisfies [6,30] 

042 = 0, 044/040 = ±\J 5/14. (5) 


The nuclear surface (1) then reduces to 


R(0,<f>) = R o 


^1 + 040 


>4o(M)± 


JUY4-4 (M) + n +4 (A 



( 6 ) 


Such shapes cannot be generated with the restricted parameterization (3), 
which therefore is insufficient for the present purpose. For positive values of 
0.40 an octahedron is obtained while for negative a 40 one finds a cube, the dual 
of the octahedron, both shapes having octahedral symmetry. The sign of the 
ratio 044/040 does not affect the intrinsic shape; the opposite sign corresponds 
to the same shape rotated over 7 t /2 around the z axis. The four different cases 
are illustrated in Fig. 1. 

In the parameterization (4) the first of the octahedral conditions (5) implies 
64 = 0, 84 = 7T or 74 = 0. For 84 = 0 or 84 = 7 r, the second of the octahedral 
conditions (5) is automatically satisfied for any value of 74 , resulting in a 
positive ratio 044 / 040 ; if 84 = 0 , the shape is an octahedron while for 84 = 
7 T it is a cube. For the remaining case of 74 = 0 one finds an additional 
solution with octahedral symmetry for $4 = arccos(l/ 6 ), corresponding to a 
rotated octahedron, see Fig. lc, but no rotated cube can be obtained with the 
parameterization (4). 

In the quadrupole parameterization (2) a given couple ((3 2 , 72 ) with /3 2 € 
[0, Too [ and 72 G [0, 7 t/ 3 ] corresponds to a unique intrinsic shape. This is not 
the case for the hexadecapole parameterization (4), that is, different triplets 
(/? 4 , $ 4 , 74 ) with At £ [0,Too[, A 4 G [0,7r] and 74 G [ 0 , 7 t/ 3 ] may lead to the 
same intrinsic shape, differently oriented with respect to the laboratory frame. 
An example of the latter are the triplets (An&i — 0,74 = anything) and 
(/? 4 , A 4 = arccos(l/ 6 ), 74 = 0 ) which correspond to the same intrinsic shape 
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(a) octahedron 


(b) cube 




(c) octahedron 


(d) cube 



Fig. 1. Hexadecapole deformed surfaces with octahedral symmetry. Four cases are 
shown: (a) and (c) are octahedrons with 040 = 0.2; (b) and (d) are cubes with 
040 = —0.2. Octahedral symmetry is obtained for 042 = 0 and 044/040 = -F-^/5/14 
[(a) and (b)], or 042 = 0 and 044/040 = — -y/ 5 / 14 [(c) and (d)]. The two signs of 
044/040 correspond to the same intrinsic shape but rotated over 7r/2 around the z 
axis. In terms of the parameters in Eq. (4): (a) 04 = 0.2 and 5 4 = 0, (b) 04 = 0.2 
and 64 = it, (c) 04 = 0.2, 84 = arccos(l/ 6 ) and 74 = 0, (d) not possible. 

for any value of 04. Shapes that are intrinsically the same lead to identical 
conditions on the interaction parameters in the Hamiltonian, as will be shown 
below. 


3 The sdg interacting boson model 


Since the bosons of the IBM represent pairs of valence nucleons, a nucleus 
is characterized by a constant total number of bosons N which equals half 
the number of valence nucleons. An important feature of the sd-IBM is the 
existence of a U( 6 ) dynamical algebra, the substructure of which leads to 
analytically solvable limits, also called dynamical symmetries. The sd-IBM can 
successfully describe quadrupole collective states in even-even nuclei but other 
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features require an extension of the sd-IBM. In particular, the hexadecapole 
degree of freedom requires the introduction of a g boson (£ = 4) and the 
upgrade of the dynamical algebra from U(6) to U(15). We do not cite here 
the many papers related to the sdg -IBM but refer the reader to the excellent 
review by Devi and Kota [31]. 


3.1 Hamiltonian of the sdg-IBM 


Since the Hamiltonian of the sdg-IBM conserves the total number of bosons, 
it can be written in terms of the 225 operators b\ m bg' m ' where b\ m (be m ) creates 
(annihilates) a boson with angular momentum £ and z projection m. This set 
of 225 operators generates the Lie algebra U(15). A boson-number-conserving 
Hamiltonian with up to two-body interactions is of the form 

H = Hi + H 2 . (7) 


The first term is the one-body part 


H\ = ejs 1 " x §] (0) + e d V / 5[d t x d] (0) + e g \ /^[g 1 * x g] (0) 

= e s s' • s + ed d' ■ d + e g g^ ■ g 

= e s n s + e d h d + e g h g , ( 8 ) 

where the multiplication x refers to coupling in angular momentum (shown 
as an upperscript in round brackets), the dot • indicates a scalar product and 
b tm = H'-V™. Furthermore, fq is the number operator for the i boson 
and the coefficient q is its energy. The second term in the Hamiltonian (7) 
represents the two-body interaction 


Ho = 


) Lv e 1 e 2 e ' 1 e , 2 




\Jif- + 


X 4] (i) • Ih X b,] {L \ 


(9) 


where the coefficients ^ are the interaction matrix elements between 
normalized two-boson states, 


v 


L 


(£i1 2 -LM l \H 2 \£'/ 2 -LM l ). 


( 10 ) 


It is henceforth assumed that t\ < i 2 and l\ < 1' 2 . 

Once the single-boson energies and interaction matrix elements vf ^ are 
known, the most general two-body sdg-Hamiltonian is uniquely determined. 
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3.2 Classical limit of the sdg-IBM 


The classical limit of any boson Hamiltonian is defined as its expectation 
value in a coherent state [32], This yields a function of the deformation vari¬ 
ables which can be interpreted as a total energy surface depending on these 
variables. The method was first proposed for the sd-IBM [33,34], The exten¬ 
sion to the sdg-IBM was carried out by Devi and Kota [35] for the simplified 
parameterization (3). 

The generic form of the coherent state for the sdg-IBM is 

I N; a 2jU , a 4#i ) oc T(a 2jU , o), ( 11 ) 


where 


r(a 2M , a iff) — s^ + ^2 a 2 ^i (d-n + d+ M ) 

n= 0,2 


+ X! a 4n(g- ll + g+ti), 

11=0,2,4 


( 12 ) 


and |o) is the boson vacuum. The a a m have the interpretation of shape variables 
appearing in the expansion (1). Since the deformation in the IBM is generated 
by the valence nucleons only, in contrast to the geometric model of Bohr and 
Mottelson [1] where it is associated with the entire nucleus, the shape variables 
in both models are proportional but not identical [36]. In terms of the variables 
in the parameterizations (2), (3) and (4), fd 2 and /J 4 are proportional in both 
models while the angles 7 or y 2 , 74 and £4 have an identical interpretation. 

The coherent state based on the parameterization (3) 


IV A, A, 7) = 


r(ft,A,7)», 


JV!(i + A 2 + A 2 ) N 


(13) 


with 


T(/3 2 ,/3 4 ,7) = s t + fd 2 [cos 7 dp + ^/§sin 7 (dl 2 + d\_ 2 )\ 

+/ 3 4 \ l (5 cos 2 7 + 1 )gl + yif sin 27 (gl 2 + g\_ 2 ) 


+ V# sin 7 (9 -4 + 9+a ) 


(14) 


was used by Devi and Kota [35] to derive the geometry of the different 
limits of the sdg-IBM. The expression for the classical limit of the general 
sdg-Hamiltonian (7) with this coherent state was given in Ref. [37]. 



If one is interested in octahedral shapes and how they appear in the sdg-IBM, 
the general parameterization (4) is needed, and the appropriate coherent state 

is 


l-W; P2, At, 72, 74 , £4) 


r(^2,^4,72,74,^4) iV |o), 


N\{l + ft + pl) N 


(15) 


with 


r(/3 2 ,/?4, 72 , 74 ,^ 4 ) = s t + /?2 cos 72 do + ^1 sin 72 ( 1 / 2 + d\ 2 ) 
+/?4 ( \[^2 COS ^4 + \[^2 sin A COS 74)^0 
- sin 5 a sin 74(^12 + 9 + 2 ) 


(16) 


+ f \!h cos h 4 - \/J sin ^4 cos 74) (gl _ 4 + g\_ 4 ) 


The classical limit of any Hamiltonian of the sdg-IBM is, for the general 
coherent state (15), defined as 

(H) = (N ;&, An 72 ,74, 6 4 \H\N;p 2 ,p 4 , 72 ,74 A}- (17) 


Once the form of the coherent state has been determined, the expectation 
value (17) can be obtained with the method of differentiation [38] which lends 
itself ideally to programming in a symbolic language. The classical limit of the 
one-body part (8) is 


£s + + e s /3f 

1 + ft + Pl 


(18) 


while that of its two-body part (9) can be written in the generic form 




N(N-l) 

W+W+W 




c ki + E c ki cos(i7 2 + A) AwA) , 

ij 


(19) 


where the coefficients Cki and c'h can be expressed in terms of the interactions 
m ■ The expressions for the non-zero coefficients Cki are 


COO 2 v ssssi 

C40 = ^dddd + 
r , _ 1 .,0 
C22 3^5 V ddgg 


C20 — 

l v 2 
7 u dddd 


1 7 ,° 

5 u ss-dd 


+ V 


sdsdi 


c 02 — 


_i_ — 

^ 35 u ddddi 


+ 


_7 4 1 1„.2 1 _4_,4 

^715 dd-gg ~ 6 dgdg ~ \\ L dgdg 


1 

3 


2° 

ss-gg 


+ V 


4 

sgsgi 


+ —V 
' 6 U 


5 

dgdg 


, 10 6 
' 33 L dgdg ) 
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( 20 ) 


„ J_ 0 , 38 2 , _89_ 4 , _62_ 6 , 1129 8 

04 18 L 9999 ' 693 L 9999 ' 1001 L 9999 ' 495 l 3S33 ' 6435 1 3333’ 

while those for the non-zero coefficients are 


„30 _ 
c 30 — 


„oo _ 
C 21 — 


_J^J 2 
y ^7 u sd-ddh 


1 7 ,2 4 - 
3 V sd-dg ' 


— 7 ? 4 

5 17 sg-ddh 


21 /10 2 , /3 4 

C 21 — V 21 U sd-dg + \j -jUgg-ddi 


3„,4 


r 1 " 1 

c 12 


= —9)2 

3 V 11 sd-gp 


4 4 


r 12 = 2 /X 7 ,2 I 4 7 ,4 
12 3 V 77 sd-gg ~ ^385 s 9'dgi 


„00 

c 03 — 


_ 

V429 *9-99 


„03 _ 
c 03 — 


4 1 

„00 _ 
l 22 — 


7 V 3 V dd dg 7\/ll V dd-dgi 


/ 10 ,- 
3 


^ 4 

3003 sg-ggi 

15 „,.4 


«30 _ 


^31 


2 „.2 


= ~\ tttV 


21 ^ dd-dg 


- „ /f 2 ? , 4 

V 77 17 dd-dg 1 


v / 7l5 l ' dd -gg 
'SL ,2 


r 2-2 

c 22 


r 21 - 
c 22 ~ 


-4 -hyO 

it; ^ 


42 V dgdg ~7~ 15 u dgdg 385 U dgdg " 1_ 30 U dgdg 


_27 4 , J_ 5_L_ 

Qfit; u dndn 1 r>r\ u dndn 


J_7,6 

165 L dgdg’ 


2_ 

'21 


11 v dd-gg 7 


2 /10 2 , 
3 V 77 ^cW-33 “ l ” 



1 a/35 ,.,6 
+ 33 Vd 


els 1 


dgdg ’ 


= -A / 1„2 

21V 33 L dg-gg 


J _ 7 )3 I ffi-i ; 4 

15 dgdg ^11 L dgdg 


_ + 1 /5 2 _8J4 

1001 c dd-gg ^6 V 7 dgdg Us/35 dgdg 


_ I? /X,, 4 _ A /10 ? ,6 

77 V 13 L dg-gg 33V 3 L d3-33’ 


165 dgdgi 


7 7.5 

V dgdg 


r 12 - 
l 13 — 

„00 _ 
l 04 — 

„03 

l 04 — 


_ 4 _ i,2 _ 12 

3 V231 dg-gg 44 

4 „,4 


^ r 4 _4 /iS 7 ,C 

455 c dg-gg 33 V 3 ^S'SS’ 


__2_ 2 _|-4_.y4 -|- -Z_<y6 

693 3333 ~ 3003 3333 1 495 3333 


/U 6 

3 1 
16 8 
6435 3333’ 


2 

‘99 


j-h_ 

7 v 9999 ' 429 


hi 4 + 

7 3333 ^ 99 


1^6 
5 3333 


16 

1287 


7y8 

5 3333’ 


( 21 ) 


where the notation i^ 2 .^ = vf li2e , e , + is introduced for (X^) 7 ^ 

(£^ 2 ) since this combination consistently occurs due to the hermiticity of the 
Hamiltonian. Furthermore, the 4)^(8 4 ) are functions defined as follows: 


030 (04 ) - 

= 1 , 0 m 

($ 4 ) = cos 

^4, 021 (04) — 

sin 84 

1 


012 1 (^4) 

= sin 2 84 

, 0i2 (^ 4 ) 

= 1 — cos 2 84 , 




0 ») = 

= 6 cos 84 

+ COS 3(54, 

003 (0i) = 3 sin5 4 - 

- sin 

3^4, 

0m 1 (^ 4 ) 

= sin 84 , 

0!?(X) = 

= COS 84 , 




0^4) = 

- 02 2 “ 2 (^; 

) = 1-cos 25 4 , 022 (^ 4 ) 

= sin 

25 4 , 


013 1 (5 4 ) 

= sin 84 + 2sin3<5 4 , 

4 >\ 3 (^ 4 ) = cos<5 4 - 

cos 

3<5 4 , 

004 (^ 4 ) = 

= 2 cos 28 , 

! + 17cos4<5 4 , 0q 4 (5 4 ) = 

2 sin 

2 84 

— sin 4 84 


These functions can be written concisely as 

0 or 1 

<t>Z( s 4) = Z c% n <pj(n 6 A ), 

n=l,l— 2 ,... 


( 22 ) 


(23) 


where (pj(9) is cos 9 (sin 0 ) for even (odd) j, with the following interaction- 
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independent constants: 


r 30 

c 300 

= 1, 

r 00 - 
c 211 — 

1, 

r 21 - 1 

c 211 — L i 


r 1 " 1 

c 120 

= 0 , 

r 1 ” 1 - 
c 122 _ 

= 1, 

r 12 - 1 
c 120 — H 

r 12 - -1 

l 122 ~ L i 

r 00 

l 031 

= 6, 

r 00 - 
c 033 — 

1, 

r 03 - 3 

c 031 — 

r 03 - -1 
c 033 — L i 

r 1 " 1 

l 311 

= 1, 

r 30 - 
c 311 — 

1, 



„oo 

L 220 

- c 2 " 2 

= 1, 

„oo 

l 222 

- c 2 - 2 - 

~ 1) c 220 = c 222 = 

c 1 " 1 

l 131 

= 1, 

C 133 = 

= 2, 

c 12 - 1 

L 131 — H 

c 12 - -1 

l 133 — H 

J)0 

l 040 

= 0, 

„00 ... 
l 042 — 

2, 

r 00 — 17 

l 044 — 1 '? 

«03 n 03 _ 9 03 

l 040 — u > l 042 — l 044 


-1. (24) 

The classical limit of the total Hamiltonian (7) can therefore be written as 


(H) = E(0 2 ,fa, 72,74,^) 


N(N-l) 

(1 + (5 2 + W 


£02 04 


c'kl + £ c ki cos(f 7 2 + J 7 4 ) 0 w( 5 4 ) , 
ij 


(25) 


where c' kl are the modified coefficients 


Coo — c oo + c' s , 

c 40 = c 40 + Q, 


c 2 o — c 2 o + <4 + e ' d , c ' 02 — c 02 + e' s + e' g , 

c 22 = C 22 + Crf + Cg, Cg 4 = C 04 + e ' g , 


(26) 


in terms of the scaled boson energies e' e = e^/(IV — 1). While the differentiation 
technique [38] allows a secure derivation of the expectation value (17), the 
particular representation (25) in terms of functions cos(i 7 2 + jli) 
is not obtained automatically. The correctness of the latter representation can 
be proven by use of trigonometric conversion algorithms which show it to be 
identical to the expression found with the brute-force differentiation technique. 

The quantum-mechanical Hamiltonian (7), if it is hermitian, depends on three 
single-boson energies q and 32 two-body interactions . In the classical 

limit with the most general coherent state (15), the number of independent 
parameters in the energy surface E((3 2 -,f3 4 , 72 , 74 ,^ 4 ) is reduced to 22 (six co¬ 
efficients c' kl and 16 coefficients c kl ). For comparison, if the simpler coherent 
state (13) is taken [37], this number is further reduced to 15. 
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4 Octahedral shapes in the sdg-IBM 


4-1 Principle of the method 


It is clear that a catastrophe analysis of the energy surface E (/3 2, Pa, 72,74, £4) 
with its five order parameters and 22 control parameters is beyond the scope 
of any reasonable analysis. Fortunately, this is not needed if one is interested 
in the realization of octahedral symmetry in the sdg-IBM. For this purpose 
one just wants to know what are the conditions on the interactions in the 
•S'dry-Hamiltonian for the surface (25) to have a minimum with octahedral 
symmetry. As shown in Sect. 2, a shape with such symmetry occurs for (i) 
P2 = 0, Pa fz 0, 72 = anything, 74 = anything and £4 = 0, (ii) p 2 = 0, 
P 4 7^ 0, 72 = anything, 74 = anything and 64 = 7r, or (iii) p 2 = 0 , /3 4 fz 0, 
72 = anything, 74 = 0 and 5 4 = arccos(l/6). 

The conditions for the energy surface E(j 3 2, /id, 72,74, 64) to have an extremum 
at p* are 


d E 

dE 

dE 

dE 

dE 

d P 2 

p * ~ dp A 

p * d 3/2 

p* ~ ^4 

* dd 4 

p* ^ 


where p* = ((3PI, 72 , 74 , £ 4 ) is a short-hand notation for an arbitrary critical 
point. Furthermore, a critical point with octahedral symmetry shall be denoted 
as o* which implies that o* is one of the three cases (i), (ii) or (iii) listed above. 
While the Eqs. (27) are necessary for E(fi ’ 2 , An 72 , 74 , &i) to have an extremum 
at p*, the conditions for a minimum require in addition that the eigenvalues of 
the stability matrix [he., the partial derivatives of E(/3 2, [3 4 , 72,74, £4) of second 
order] arc all non-negative. All expressions in this section are obtained starting 
from the generic expression (25) for the energy surface and its derivatives up 
to second order. 


4-2 Extrema with octahedral symmetry 


Let us now apply the above procedure to the case of octahedral symme¬ 
try which requires the establishment of an extremum of the energy surface 
E{@ 2 , Pa, 72,74, fa) for P* = o*■ 

Consider first the case of an octahedral shape (i) or a cubic shape (ii), cases 
that can be treated simultaneously. Four of the five extremum conditions (27) 
are identically satisfied for p* = o* and do not lead to any constraints on the 
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coefficients c' kl and c kl . The derivative in /3 4 leads to the equation 

PI [-44 + 2 cq 2 ± 214/31 - (2 c' 02 - 44 - 76 cq 4 ) / 3 f T 7 c“/ 3 f] = 0, ( 28 ) 

where the upper (lower) sign applies to <5 4 = 0 (<5 4 = 7 r). The spherical point 
/3 4 = 0 is always an extremum of the energy surface. Other extrema /3 4 are 
found as solutions of a cubic equation and therefore have cumbersome expres¬ 
sions. In the search for octahedral minima one is not interested in numerical 
values or analytic expressions for /3 4 but one simply wants to know whether a 
hexadecapolc deformed minimum exists or not. This question can be readily 
answered if one assumes the coefficients of the odd powers of /3| in Eq. (28) 
to be zero, which happens for 4 = 0. In that case, the non-zero solutions of 
Eq. (28) are 



leading to the conclusion that a real, positive /3 4 is found if the combinations 
24 — Cq 2 and —4 + 2c ' 04 + 384 have the same sign. For non-zero values 
of CQ 3 [which is related to the s-g mixing matrix element vj a . gg , see Eq. (21)] 
the analysis is more complicated. The cubic equation (28) has real coefficients 
and therefore it always has at least one real solution /3|. In addition, since the 
ratio =FCo 3 /(—4cg 0 + 24 ) is negative for one of the choices <5 4 = 0 or 5 4 = n, 
it follows that the solution /3 4 is positive in that case. We conclude that there 
exists always an extremum with octahedral symmetry for any sdg -Hami 11 onian 
except in the pathological case of no s-g mixing, xP = 0 , in which case the 
condition is that the combinations 2 cq 0 — c ' 02 and —c ' 02 + 2 c ' 04 + 384 should 
have the same sign. 

For the case (iii) with <5 4 = arccos(l/ 6 ), the derivatives in y 2 and q 4 are 
identically zero and do not lead to any condition. The derivatives in /3 2 , /3 4 
and 5 4 lead to the equations 

Pt 2 [-3 (Vttclz 1 + 354) + 7 (v^cij 1 - 54) /?*] = 0 , 

PI [324 (-24 + 4) + 63 (4c“ + 5 a/354) PI 
+4 (-814 + 162cq 4 + 18534 + 35^4) Pf 
- 21 (44 + 5V354) pf] = 0, 

pf [9 (- 2 V 354 + 7 4) + 224 (^4 - 4) Pi] = 0. (30) 

It would therefore seem that for the energy surface (25) the critical conditions 
in the case (iii) lead to equations that are different from the those obtained 
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in the cases (i) and (ii). It should not be forgotten, however, that the coef¬ 
ficients c' kl and c kl are expressed in terms of single-boson energies q and in¬ 
teraction matrix elements e . After substitution of the expressions given 
in Eqs. (20), (21) and (26), it is found that the Erst and third conditions of 
Eq. (30) are identically satisfied while the second reduces to the one given in 
Eq. (28) with the upper sign. This confirms the earlier statement that a given 
intrinsic shape leads to unique conditions on the single-boson energies and 
interaction matrix elements, independent of the orientation of that shape in 
the laboratory frame The result provides an additional and independent check 
on the correctness of all the equations involved in this comparison. Therefore, 
the analysis henceforth can be restricted to the cases (i) and (ii) of octahedral 
and cubic intrinsic shape. 


4-3 Minima with octahedral symmetry 


So far, no constraints are found on the single-boson energies q and interaction 
matrix elements vf e e i g i since Eq. (28) has always a real, positive solution, 
either for h 4 = 0 or for o* 4 = 7 r, except in the pathological case of no s-g mixing 
mentioned in Sect. 4.2. Constraints are found by requiring that the extremum 
is a minimum or, equivalently, that the eigenvalues of the stability matrix are 
all non-negative. No conditions follow from second derivatives involving y 2 and 
74 at a critical point o* with octahedral symmetry. Furthermore, the second 
derivatives involving (3 4 are decoupled from those pertaining to /3 2 and <5 4 , that 
is, the following equations are identically satisfied: 


d 2 E _ d 2 E 
dfcdfa o . d(3 4 dS 4 


(31) 


so that the stability in (3 4 follows from the inequality 


d 2 E 

M 


> o. 

o* 


(32) 


Some insight can be obtained by assuming the odd powers of (31 to be zero, 
eg® = 0, in which case the condition (32) reduces to 

(2 C 00 — c 02)( c 00 — c 02 + c 04 + 19 c 04) ^ , 

-c ' 02 + 2cq 4 + 38 c™ - U - 


Since the combinations 2c ' 00 — c ' 02 and —c ' 02 + 2c ' 04 + 38c[] 4 must have the same 
sign (see Sect. 4.2), it follows that the conditions 

2cqo — c 02 — 0) ~ c o2 A 2c 04 A 38cq 4 > 0. (34) 
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are necessary and sufficient for the energy surface E(f3 2 ,fi 4 , 72 , 74 ,^ 4 ) to have 
an extremum at non-zero hexadecapole deformation which is stable in 6 a- 

The stability in 62 and <5 4 follows from the diagonal derivatives, 


d 2 E 


d 2 E 

W 


4cqo + 2c2q 4= 2c2?/3| (4cq 2 2c2 0 < ^c' 22 )6a 2 


T(28c°° - 2cf 1 )6f - (4c' 04 - 2J 22 + 76cg°)^ 4 
(1 + /9I 2 ) 3 

Tl5 c° 0 l/3f - 280 c™6t 

(i + 6f) 2 


(35) 


These equations are coupled, however, since the off-diagonal derivative gener¬ 
ally is non-zero, 


d 2 E 

d[3 2 dS 4 


2ci2 1 /^4~ 4= 7c\ 3 1 6l 

(1 + 6 f ) 2 


(36) 


Again, simplifications arise if the coefficients of the odd powers of 64 vanish, 
Cq° = C 21 = 0 . If, in addition, the off-diagonal derivative vanishes, cj ^ 1 = 
C 43 1 = 0 , the stability in 62 and 84 at the hexadecapole deformation 61 given 
by Eq. (29) is guaranteed by the following conditions: 

( c 02 — c 2 o) ( 1c 02 + c 22 ) + (2c 00 — c 20 )(2c 04 — c 22 + 38c 04 ) < 0, (37) 


and 


c 


00 

04 


< 0 , 


(38) 


respectively, where use has been made of the conditions (34), required to have 
an extremum which is stable in 6 a- 

Provided that c °3 = c 24 = c \ 2 1 = cj ^ 1 = 0) Eqs. (34), (37) and (38) are the 
necessary and sufficient conditions for the energy surface E( 6 2 , 6 a , 72 , 74 , 84 ) to 
have a local minimum with octahedral symmetry. It is not necessarily a unique 
minimum and it may not be the global one. In particular, the conditions do 
not exclude the existence of a quadrupole-deformed minimum with 62 7 ^ 0 at 
some hexadecapole deformation which differs from 6 a given in Eq. (29). To 
exclude the latter possibility for all 64 , the following stronger conditions must 
hold [see the first of Eqs. (35)]: 

2cqo c 20 — 2c 02 — c 20 — c 22 A 0, 2c 04 — c 22 + 38c[] 4 < 0, (39) 
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Fig. 2 . Illustration of the energy surface ^4,72,74, £4) of Eq. ( 25 ) with co¬ 

efficients Cqq = — 0 . 2 , C 20 = — 0 . 5 , Cg 2 = — 1 , C40 = 0 , c '22 = — 0 . 25 , Cg 4 = 0.15 
and Cq4 = — 0 . 02 , in arbitrary energy units times N(N — 1 ). Left: (/ 32 ,A 0 plot for 
72 = 74 = <$4 = 0 . Right: (,84,54) plot for = 0 and 72 = 74 = 0 . Blue areas 
correspond to low energies. 

still assuming that C 03 = = 0. These stronger conditions can be combined 

with the inequalities (34) to yield, with the help of Eq. (26), 

c 02 + e g A 2cqo < C 20 + AO C 02 < 2 cq4 + 38cq4 + e” < C 22 + AO (40) 


where e" are the scaled single-boson energies relative to the s boson, e" = 
e'z — e' s . These conditions have the advantage to be sufficiently simple to allow 
an intuitive understanding. With reference to Eqs. (20) and (21), the inequal¬ 
ities (40) express the condition that (i) the s-g mixing is strong enough as 
compared to the energy difference e g — e s to develop a hexadecapole deformed 
minimum and (ii) the s-d and d-g mixing is sufficiently weak as compared to 
the energy differences e s — ej and e g — e^, respectively, so as the minimum to 
remain at 62 — 0 . 

I 11 Fig. 2 an example is shown of an energy surface E((3 2 ,f3 4,72,74,^4) with 
coefficients in Eq. (25) that satisfy the conditions (34), (37) and (38). It is 
obviously not possible to display a five-dimensional surface in its full com¬ 
plexity and judiciously chosen two-dimensional intersections must be shown 
to illustrate the structure of E(/3 2 , 72, 74, $4). In the study of shapes with 

octahedral symmetry one may take 72 = 74 = 0 since there is no dependence 
on these variables at the minimum. Figure 2 shows two intersections, in the 
(^ 2 , AO plane with <f 4 = 0 and in the (A 4 , <f 4 ) plane with /?2 = 0 , and confirms 
the existence of two minima with octahedral symmetry for Aj = 0 and ^ 0 , 
with 51 = 0 and 5 4 = ir, corresponding to an octahedron and a cube, respec¬ 
tively. A third minimum is seen at 5 4 = n/2 which represents a shape with a 
lower discrete symmetry. 
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5 Conclusions 


In this paper a general Hamiltonian of the sdg-IBM with up to two-body 
interactions was analyzed as regards the question of the occurrence of intrinsic 
shapes with octahedral symmetry. Such an analysis requires the use of the 
most general sdg coherent state (16), leading to a classical energy surface of 
the generic form (25). A stability analysis of this surface leads to a set of 
conditions, Eqs. (34), (37) and (38), which are necessary and sufficient for the 
occurrence of a minimum with an intrinsic shape with octahedral symmetry. 

Due to the complicated nature of the stability conditions, only qualitative 
conclusions have been drawn at this point with regard to the occurrence of 
octahedral shapes in the sdg-IBM. More concrete conclusions will be drawn in 
the study of a simpler .srig-Hamiltonian, which has the advantage of exhibiting 
dynamical symmetries and which is the topic of another paper in this series. 
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